lammps教程：python ovito编程绘制速度分布 |
您所在的位置：网站首页 › lammps数据处理 python › lammps教程：python ovito编程绘制速度分布 |
lammps教程：python ovito编程绘制速度分布
|
大家好，我是小马老师。 本文介绍使用python和ovito绘制速度分布曲线的方法。 在lammps流体模拟中，都会用到速度分布曲线，以应流体的流动性能。 有两种方法绘制速度分布曲线，一是在lammps in文件中使用分块计算的原理直接输出速度分布数据，二是使用编程的方法对数据进行后处理得到速度分布曲线。 本文介绍第二种方法，使用python和ovito绘制速度分布曲线。 ovito处理速度分布和lammps原理基本一致，也是先进行分块切割，得到每一块的原子数据后，再进行数据求和平均。 使用的命令为：SpatialBinningModifier 这个模块在ovito软件中也有，不过是pro版本才可以用。 如果没有ovito pro版本，则可以使用python调用ovito库函数的方法实现同样的功能。 这种方法完全免费，只需要安装python和ovito库即可。 安装方法： python ovito模块安装方法 安装完成后，新建flow.py，复制粘贴下面的代码。 该案例模拟了原子沿x方向的流动过程，需要计算的是原子沿y方向的速度分布。 from ovito.io import import_file, export_file from ovito.modifiers import SpatialBinningModifier #读入dump文件 pipeline = import_file("dump.xyz") #分块处理速度数据 pipeline.modifiers.append(SpatialBinningModifier( property = 'Velocity.X', direction = SpatialBinningModifier.Direction.Y, bin_count = 15, reduction_operation = SpatialBinningModifier.Operation.Mean, )) #输出 export_file(pipeline, 'velocity_profile.txt', 'txt/table', key='binning')运行python： python flow.py运行完成后，会生成一个文本文件：velocity_profile.txt，该文件一共两列数据，第一列是Y方向的每个分块bin的坐标，第二列为原子X方向的速度。 绘制曲线如下图所示： |
------------------------------- 公众号：lammps加油站【本文地址】
今日新闻 |
推荐新闻 |